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We study the zero-temperature magnetization curve (M — 
H curve) of the one- dimensional quantum antiferromagnet of 
spin one. The Hamiltonian H we consider is of the bilinear- 
biquadratic form: H = /(sVsi+i) (+Zeeman term) where 
Si is the spin operator at site i and f{X) — X + f3X 2 with 
< /3 < 1. We focus on validity of the 5- function bose- 
gas picture near the two critical fields: upper critical field 
H s above which the magnetization saturates and the lower 
critical field H c associated with the Haldane gap. 

As for the behavior near H a , we take "low-energy effective 
S-matrix" approach where correct effective bose-gas coupling 
constant c is extracted from the two-down-spin S-matrix in 
its low-energy limit. We find that the resulting value of c 
differs from the spin- wave value. We draw the M — H curve 
by using the resultant bose gas, and compare it with numeri- 
cal calculation where the product-wavefunction renormaliza- 
tion group (PWFRG) method, a variant of the S. R. White's 
density-matrix renormalization group method, is employed. 
Excellent agreement is seen between the PWFRG calculation 
and the correctly-mapped bose gas calculation. 

We also test the validity of the bose-gas picture near the 
lower critical field H c . Comparing the PWFRG-calculated 
M — H curves with the bose-gas prediction, we find that there 
are two distinct regions, I and II, of j3 separated by a critical 
value f3 c (~ 0.41). In the region I, < ft < f3 c , the effective 
bose coupling c is positive but rather small. The small value 
of c makes the "critical region" of the square-root behavior 
M ~ y/H — H c very narrow. Further, we find that in the (3 — > 
/3 C — 0, the square-root behavior transmutes to a different one, 
M ~ (H-H c ) e with 9 « 1/4. In the region II, /3 C < < 1, the 
square-root behavior is more pronounced as compared with 
the region I, but the effective coupling c becomes negative. 

PACS numbers: 75.10.Jm, 75.40.Cx, 75.45.+j 



I. INTRODUCTION 

The magnetization process (M — H curve, M: mag- 
netization, H: magnetic field) of one-dimensional (ID) 
quantum spin system has recently drawn much atten- 
tion, due to the remarkable progress in material synthe- 
sis techniques and high-field experiments. Q Spin chains 
with various S (spin magnitude) and/or with non-trivial 
special structures, etc, exhibit various interesting mag- 
netic behaviors (e.g., field- induced phase transitions, like 
plateau in the M — H curve 0), many of which still await 
theoretical analyses. 



As for the M — H curve of the gapful S = 1 antifer- 
romagnetic (AF) Heisenberg chain, it has been known 
that, on raising the magnetic field from zero, there is a 
critical field H c above which the system becomes magne- 
tized. This critical field relates to the excitation gap (Hal- 
dane gap) A as H c — A/(g/iB) (, g: g-factor, /ib: Bohr 
magneton) . We call H c lower critical field because there 
also exists upper critical field H s (saturation field) above 
which the magnetization saturates to M s = 1. Near H s , 
it is well established that the M — H curve behaves as 
M s - M - y/H s - H. JT|-[l4[] Near H c , similar behavior 
M ~ (H - H c f with 6 « 1/2 has also been known, ffl 
but numerically, whether the exponent 8 is exactly 1/2 
or not, has remained to be less conclusive as compared 
with the behavior near H s . 

The expected square-root behavior M ~ \JH — H c 
has been explained via approximate mapping to the <5- 
function bose gas |7|§] or to the fermion ga s. There 
is numerical evidence for such mappings, [|s|,[l0[ but, the 
square-root behavior of the bulk magnetization itself has 
not been fully verified yet. In the finite-iV (N: sys- 
tem size) diagonalization study, || the smallness of N 
disables us to make quantitative discussion of the bulk 
magnetization near H c . In Ref. Jic|| , a large-size system 
is treated by the density-matrix renormalization group 
H (DMRG, for short) to study the low-lying excitations 
having small 5 t z ot (< 4). It was shown there that the low- 
lying excitations admit fermionic interpretation, just as 
in the case of the Bcthe-ansatz solution of the 5-function 
bose gas at low particle density. This result is indeed a 
strong support for the bose-gas picture, but the smallness 
of St ot implies that, in a strict sense, the result applies 
only to the system in the vanishing magnetization density 
m = SfjN -> 0. 

Further, there is a quantity which is undoubtedly im- 
portant for the bose-gas picture, but has not been con- 
sidered seriously so far: the effective bose-gas coupling 
constant c. Let us consider a ID (effective) (5-function 
bose gas with the Hamiltonian 

/* w t ( , W( , )+e#W( , W( , W( , )1 , 

(i.i) 

with (x) and (j>(x) being the bose operators. By solving 
the Bethe ansatz integral equation we obtain the ground- 
state energy density e as [O 
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r = p: particle number density. 
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In terms of the particle density p, we have 
e(p) = A 3/9 3 + A 4 p 4 + • • • , 

with 

A 3 = CT7T 2 /3 

A 4 = -4(rV/(3c). 



(1.3) 



(1.4) 



(1.5) 



(1.6) 



In the bose-gas picture, the M—H curve of the spin chain 
near the critical field corresponds to the p — p curve of 
the bose gas, where p = de(p) /dp is the chemical poten- 
tial. As far as the critical behavior is concerned, only 
the p 3 -term is relevant, which determines both the crit- 
ical exponent (= 1/2) and the critical amplitude. Since 
the coefficient A3 does not depend on c, the critical be- 
havior must be c-independent and "universal". However, 
actual p — p curve heavily dependents on the value of c 
which comes from p 4 and/or higher-order terms. In fact, 
the coefficient A\- of p fc -term (k > 4) is proportional to 
c 3 ~ k , which becomes dominant for small c. Accordingly, 
the critical region of the square-root behavior will be- 
come rather narrow for small c. Hence, knowledge of the 
actual value of the effective coupling constant c is indis- 
pensable, in order to make a fully quantitative test of the 
bose-gas picture. 

The aim of the present article is to give a quan- 
titative test for the bose-gas picture of the bilinear- 
biquadratic AF chain in a field. For this purpose, we em- 
ploy the quantum version |3j of the product- wavefunction 
renormalization group method Q (PWFRG method, for 
short). The PWFRG is a variant of the S.R. White's 
DMRG, H which is specially designed to obtain the 
"fixed point" (=thermodynamic limit of the system |J) 
of the DMRG iterations efficiently. In Ref. §, it is shown 
that the PWFRG which was originally implemented for 
2D classical systems, can also be applied to ID quan- 
tum systems by replacing the transfer-matrix multipli- 
cation with the modified Lanczos operation. Even with 
relatively small number of retained bases, which is con- 
ventionally denoted as "m" , PWFRG calculations accu- 
rately reproduce exact M — H curves for integrable mod- 
els, |^Jl^l which demonstrates both the efficiency and the 
reliability of the method. 

The hamiltonian of the model we consider is 



i-zccman; 



^BLBQ = ^ \S*i ■ Si + ± + 0(Si ■ Si + i) 2 



where gps = 1, or, these factors are absorbed into the 
field if. In this paper we consider the case < j3 < 1. 
We concentrate on the behavior of the M — H curve near 
the two critical fields if, and H r . 



II. NEAR THE SATURATION FIELD H s 

Let us consider first the system near the upper critical 
field (saturation field) H s = 4, where the M — H curve 
shows the square-root behavior: 1 — M ~ \/H s — H . As 
compared with the one near H c , this behavior itself is 
well-established (by exact diagonalization JlJ] and Bethe 
ansatz |l^,|l^]). Our concern is how well the (5-function 
bose gas can describe the M — H curve away from the sat- 
uration field. For this purpose, we should derive the cor- 
rect effective bose gas hamiltonian, which we shall make 
by employing the low-energy effective S'-matrix approach. 

Above the saturation field H s , the system is ferromag- 
netically ordered ("all-up" state). Low-energy excita- 
tions slightly below H s are well described in terms of 
"down spins" in the sea of up ("+1") spins. Regarding a 
down spin as a particle, we consider two-body scattering 



problem to obtain the exact two-body S'-matrix. [|13|J14 
The point is that in the low-energy limit the S'-matrix 
reduces, in most cases, to that of the (5-function bose gas 
with an effective coupling constant. With this correct 
coupling constant, the bose gas will give a quantitative 
description of the system near H s . 

By {|<7i, 0"2, . . . , >} (<Tj = 0,±1), we denote the S z - 
diagonal bases of the spin chain. To solve the two-down- 
spin problem we express the eigenvector \k, k' > in terms 
of the wavefunctions il>(x, y, fc, k') and f{z, k + k') as 



\k,k' > = y^tp(x,y,k,k')\ll---101---101---> (2.1) 

^ — ' x y 

x<y 

+ ^/(2,fc + fe')|l---l(-l)l--->. (2.2) 



The two-body S'-matrix S(k, k') is introduced through 
the asymptotic (y — x — * 00) behavior of the wavefunc- 
tion: 

ip{x, y, k, k') ~ e lkx e tk ' y + S(fc, k')e ik ' x e iky . (2.3) 

The eigenvalue problem HblbqI^, k' >= E(k, k')\k, k' > 
is solved to give 



E(k,k') = e(k) + e(k') 



(2.4) 



^zeeman H ^ ^ S^ , (1 



where e(k) = — 2 + 2cosfc is the one-particle (one-down- 
spin) energy. Since the one-particle energy e(fc) takes its 
minirnji^ at k — n, we put 

k = 7T + k, k = it + K 



(2.5) 



where Si = (Sf , Sf , Sf ) is the 5=1 spin operator at the 
site i. For notational simplicity, we have took the units 



and consider the "low energy limit" k, k' — * 0. In this 
limit, the S-matrix whose explicit expression has been 
given in Ref. Ilj], becomes 



2 



S(k,k') -> S {hosc) {n 



, 3/3+1 , 

rv : ~ J 



(2.6) 



where 5( bose ' ) (fc, c) is the 5- matrix of the 8- function bose 
gas with hamiltonian 7i( bosc )(l, c) in (1.1). Explicitly, we 
have 



S^ osc Hk,c) 



IK 



(2.7) 



Therefore, in discussing the low-energy properties near 
the saturation field, the bilinear-biquadratic chain ( |l.g| ) 
is equivalent to the bose-gas with the effective coupling 
constant 



3/3+1 



(2.8) 



A remark is in order. For /3 > or /3 < — 1/3 the cou- 
pling constant (^^) takes negative value. It has been 
known that the 5-function bose gas with negative cou- 
pling constant is unphysical, because the system is un- 
stable against formation of multiparticle bound states. 
In the present case of finite-5 spin chain, however, N- 
particle bound state with large N are kincmatically for- 
bidden because more than 25 "particles" cannot exist at 
a single site. Further, in the N = 2 case, the bound state 
actually exists, but, for < j3 < 1 it is a high-energy 
mode which can also be neglected, near the saturation 
field at least. Hence, the effective bose gas with negative 
coupling constant, does have a meaning in the present 
case. 

From the bose-gas energy density e(p) in (|l.4| ) and (1.5) 
with (7 = 1 and c = —(3/3 + l)//3, we obtain the M — H 
curve through 



M 



H s -H 



1-P, 
dp 



(2.9) 
(2.10) 



At arbitrary p, we can numerically solve the Bethe- 
ansatz integral equation Jll| by converting it to a ma- 
trix equation, to obtain de(p)/dp within any required 
precision. To check the validity of the bose-gas descrip- 
tion of the M — H curve, we performed the numerical- 
renormalization-group calculations. The method we em- 
ployed is the quantum version || of the PWFRG, 
which allows us to make fixed- H calculations sweeping 
the value of H giving the M — H curve M = M(H) in 
the thermodynamic limit. Fig.l shows comparison be- 
tween the bose-gas results and the PWFRG calculations, 
where we see excellent agreements for unexpectedly wide 
range of the field H. The validity of the bose-gas picture 
in the quantitative description of the M — H curve, is 
thus verified. 

Our two-down-spin 5-matrix approach is easily ex- 
tended to general-5 bilinear-biquadratic chain with the 
hamiltonian 



H 



Si ■ Si+i + P(Si 



Si+i) 2 



(2.11) 



After a straightforward calculation very similar to the 
ones in Refs. |l3|,|l4|], we obtain the effective bose-gas cou- 
pling as 



_ 1 + /3(3 - 85 + 85 2 ) 
1 - 5 + /3(3 - 95 + 85 2 ) ' 

For the pure bilinear case (/3 = 0), we have 

c= 2/(5- 1). 



(2.12) 



(2.13) 



In Fig. 2 we compare the PWFRG-calculated curves for 
5 = 1/2, 1, 3/2, 2 to b ose-g as curves with corresponding 
values of c given by (2.13). We again see satisfactory 
agreements. 

Let us give a comment on previous studies related to 
the present one . For the spin-5 "pure" Heisenberg AF 
(/3 = in ( 2.11 )), the bose-gas description near H s has 
already been made within the conventional spin-wave- 
theoretical approach. This approach gives the value 
of the bose-gas coupling constant to be 



= 2/5, 



(2.14) 



which is differe nt fro m (2.13). The spin- wave value (2.14) 
deviates from (2.13) very much at small 5 (even the sign 
disagrees at 5 = 1/2), although both are the same i n the 
large-5 limit. Having seen that the bose-gas with ( 2.13| ) 
gives the correct M — H curve, we must say that, at 
small 5, the spin- wave approach is not reliable enough 
for quantitative studies of the AF chains, even in the 
neighborhood of the saturation field. Note that, using 
the the Dyson-Maleev transformation and taking the con- 
tinuum limit, we can formally rewrite ||] the spin-chain 
Hamiltonian ( 2.11 ) wit h 3 = into the 8- function bose- 
gas Hamiltonian with (2.14). Although this transforma- 
tion seems to be exact in the operator level, there is a 
constraint on the state space: the boson number can- 
not exceed 25 at each site. This constraint amounts to 
"kinematical interaction" between the spin waves, which 
may be t he source of the disagreement between ( 2.13 ) 
and (|2~T4|). 



III. NEAR THE LOWER CRITICAL FIELD H c 

At H = the ground state is singlet and non-magnetic. 
On raising H , system still remains to be singlet upto a 
critical field H c above which the ground-state become 
magnetized. The field-induced phase transition at H c is 
a level-crossing transition between the singlet state and 
the lowest-energy triplet state (both at H = 0), hence the 
critical field H c is, in our unit, just the excitation gap 
("Haldane gap") A. Then, in the bose-gas description 
near H c , the singlet ground state should be interpreted 
as the "vacuum" , and the triplet state with 5f ot = 1 the 
"one-particle state" . 

For 8 « 0, the "one-particle" energy dispersion aj(k) 
takes its minimum at k — n. The dispersion curve around 



3 



this minimum is often assumed to be relativistic one, 
00 



w(jfc) = y/A 2 +v 2 k 2 , 



where k = k — ir and v is called spin-wave velocity. In 
the low-energy (|fc| — > 0) limit, (|3.1[) becomes 



w (fc) = A + — P. 



(3.2) 



Unlike the case of H ss H s , this "one-particle state" can- 
not be treated exactly, because the "vacuum" itself is not 
known exactly due to the non-integrablity of the system 
(except for some special values of (3 (= ±l,oo)). Ac- 
cordingly, for general (3, it is impossible to calculate the 
exact two-body S'-matrix S(k, k') which is utilized in the 
previous section to determine the effective coupling con- 
stant c. Nevertheless, if we assume the bose-gas picture 
to be held, we can "indirectly" determine the value of 
c from the M — H curve obtained by the PWFRG. To 
see whether the obtained value of c lies in a reasonable 
range or not, serves as a partial check of the validity of 
the bose-gas picture. 

Near H = H c = A, we should relate the M — H curve 
to the bose-gas energy density e(p) as 



H-Hr 



M = p, 

Hp) 
dp • 



o), 



Then, for the square-root behavior 

M ~ ^H-H c (H -> H c 
we should expect the expansion of the form, 
H = H c + 3A 3 M 2 + 4A4M 3 + 5A 5 M 4 + 



(3.3) 



(3.4) 



(3.5) 



where we have used (1.5). Since the expressio n (p.2| ) of 
t he o ne-particle energy implies a = v 2 /(2A) in (|1.1| ) and 
(1.6), we have 



A 3 = 
A 4 = 



3 6A ' 
-4vV/(12A 2 c) 



(3.6) 
(3.7) 



From (1.4)-(1.6), it is clear that the width of critical re- 
gion essentially depends on the reduced coupling constant 
c defined by (see ([L~6|)) 



c = c/(J, 
= -4A 3 /A±. 



(3.8) 



W 



(c) 
M 



0.1/| 7 |. 



(3.10) 



Then M < W$ implies I7MI < 0.1 which may be re- 
(3-1) garded as a necessary condition for the criticality. Cor- 



respondingly, we can introduce defined by 



W 



M I ' 



(3.11) 



which represents the width of the critical region in H — 

Let us now check the bose-gas prediction of the M — H 
curve by comparing it with the PWFRG calculations. 
From the H — M curve, we determine A (=H C ), a, 7 
and S in ( |3.9| ) by the least-square fitting. In Fig. 3 we 
show the PWFRG results of the H — M curves near H c 
for ,8 = 0, 1/3. The obtained values of A are 0.410 (for 
(3 = 0) and 0.699 (for (3 = 1/3), both of which are in 
good agreement with the known values 0.4105 (for the 
former @) and 0.699 (for the latter 

To verify the relation a = a = v 2 / (2 A) we need values 
of v. For f3 — using the known value |l(|[i3 v — 2-46 we 
have a = 7.38 which should be compared with a = 7.65 
obtained from the H—M curve; the obtained value of a is 
in reasonable agreement with a. For (3 = 1/3 there seems 
to be no serious numerical evaluation of v. We therefore 
consult Ref. [^9| where a variational calculation of u{k) 
beyond the single-mode approximation [^o| (which gives 
er = 5/9 = 0.555 . . .) is made; we have 



(32645 
0.5246. 



359V6529)/117522 



(3.12) 



This value is also in reasonable agreement with a = 0.487 
obtained from the H—M curve. Hence, the H—M curves 
reproduce the "one-particle quantities" in the bose-gas 
picture. Further, the coefficients 7 and 5 are estimated 
to be 



13 = : 7 = -8.65, 5 = 32.6 
/3 = 1/3 : 7 = -3.63, 5 = 7.46. 



(3.13) 
(3.14) 



Since the negative values of 7 implies the positive ef- 
fective bose-coupling constant, our PWFRG calculation 
supports the validity of the bose-gas picture for (3 = 
0,1/3. 

We should point out that, although the bose-gas pre- 
diction for the square-root behavior seems to be valid, 
the "critical region" of the square-root behavior in the 
M — H curve is rather narrow, since the obtained values 
of 7 and S are non-negligibly large. In fact, the quantity 



W^j defined by (3.10) characterizing the width of the 



If we rewrite (p.5f) as 



critical region, is very small: 0.012 (for (3 = 0) and 0.028 

H - H c = air 2 M 2 (l + 7M + SM 2 ) (+0(M 5 )), (3.9) ( for § = 1/3)- Corresponding values of defined by 

(pd| ) are even smaller: 0.010 [fi = 0) and 3.8 x 10~ 3 
a condition for the square-root criticality is \"fM\ « 1 (f3 = 1/3). 

(and also \SM\«1). By W$ , we denote width of the One notable behavior which we found in the PWFRG 

critical region in M, which we conveniently define as calculation is that, on raising (3 from 0, the H - M 
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curve becomes fla tter and flatter, or equivalently, the 
value of a in ( |3.9| ) becomes smaller and smaller; there 
seems to be a critical value (3 C (« 0.41) at which a van- 
ishes. Accordingly, the critical behavior of the M — H 
curve at H c changes from square-root type to another 
one ~ (H - H c ) e (9 m 0.25) (Fig.4). In the bose-gas 
picture, this change of the M — H curve may be under- 
stood as the vanishing of the fc 2 -term in the expansion of 
the one-particle excitation energy u/(k). Interestingly a 
qualitative change of the static structure factor S(q) has 
been found very near [3 C . Since both of these changes 
reflect changes in the ground state and the low-energy ex- 
citation mode of the system, it is likely that they have a 
common origin. 

Above (3 C , the square- root behavior reappears. How- 
ever, the coefficient 7 becomes positive (although small), 
implying negative effective coupling constant (Fig. 5). 
The square-root behavior itself becomes manifest due to 
small I7I, but the negative coupling disables us to take 
the naive bose-gas picture in this region of f3; for justifica- 
tion of the picture, we should inspect the "bound states", 
just as we did in discussing the M — H curve near the 
upper critical field H s . Note that for systems in the Hal- 
dane phase where the orientational order (characterized 
by the string order parameter p^ ]) exists, the "particle" 
is a moving domain wall separating two regions each of 
which has complete orientational order. |^4|,[19],|l| In this 
view, total S z carried by a low-lying excitation mode is 
the "height" of the wall. Then, if the wall width is narrow 
(~ one lattice spacing) , we can adopt a similar reasoning 
as in the previous section justifying the negative coupling 
bose-gas picture: formation of stable bound states will be 
forbidden due to the kinematical constraint that the lo- 
cal wall height (~ total S z , for thin wall) cannot exceed 
S. The actual situation is, however, subtle because the 
domain-wall is somewhat fuzzy (due to the zero-spin de- 
fects Ipij ) and its width may not be narrow. In this 
view, we should say that full justification of the bose-gas 
picture for (3 > (3 C seems to require further study. Nev- 
ertheless, the square-root behavior M ~ \JH — H c itself 
is confirmed by our PWFRG calculation. 

IV. SUMMARY 

In this paper, we have studied the zero-temperature 
magnetization process (M — H curve) of the S = 1 
isotropic antiferromagnetic spin chain with both the bi- 
linear and biquadratic forms of interactions in the range 
< (3 < 1 where (3 is the coefficient ratio between 
the biquadratic term and the bilinear term. Quantita- 
tive test for the bose-gas picture near the critical fields 
H s (saturation field) and H c (lower critical field) has 
been made with the help of the product-wavefunction 
renormalization-group (PWFRG) method which is a 
variant of S.R. White's density-matrix renormalization 
group (DMRG). 



Near H s we have derived the correct effective bose-gas 
coupling constant from the two-down-spin scattering ma- 
trix in its low-energy limit. The resulting delta- function 
bose gas yields M — H curves which are in good agree- 
ment with the PWFRG calculations. 

Near H c , the square-root behavior M ~ \JH — H c has 
been confirmed by our PWFRG calculation throughout 
the range of (3 studied. Here it should be noted a re- 
cent finite size scaling calculation by Sakai and Taka- 
hashi gave a consistent result for the (3 = case. |2j| 
We have, however, found two distinct regions of (3 sepa- 
rated by a critical value (3 C » 0.41. In the small (3 region, 
< (3 < (3 C , the effective bose-gas coupling c extracted 
from the PWFRG-calculated M — H curve is positive but 
small, making the critical region of the square-root be- 
havior rather narrow; it becomes narrower and narrower 
on approaching f3 c . At (3 C , the M — H curve seems to 
exhibit a different critical behavior M ~ (H — H c ) with 
9 ps 0.25. In the large (3 region, although the square-root 
behavior is more pronounced due to large value of |c|, 
the sign of c becomes negative, which sharply contrasts 
to the small-/? region. 

As regards the M — H curve of the bilinear-biquadratic 
Heisenberg chain, cusp-like singularities in the "middle- 
field" region have been known for integrable SU(iV) 
chains. |2q| Whether a similar behavior can also be found 
for general, non-integrable cases is an interesting prob- 
lem. Although we have concentrated on the behavior 
near the critical fields in the present paper, we have ob- 
tained a full M — H curve from H = to H = H s . In the 
large j3 region, we have actually found a clear cusp-like 
singularity very similar to the one in the SU (3) (Lai- 
Sutherland) model, |27| whose detailed account will be 
given in a separate paper. 

Finally we would like to remark that the bose-gas de- 
scription which we investigated in the present paper may 
not be the only one for "quantitative" description of the 
M — H curve of the AF spin chain. For example, in a re- 
cent paper, Yamamoto |2q] gave a different picture for the 
ground-state properties of the bilinear-biquadratic chain. 
Such an analysis may be helpful for clarifying nature of 
the system in the region (3 > (3 C . Also, "quantifying" 
other low-energy effective theories is an interesting and 
important problem. For this purpose, the approach we 
have taken in section III where microscopic quantities of 
the effective theory are extracted from bulk quantities 
calculated by a reliable method, like the DMRG. 
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FIG. 1. Comparisons of the effective (5-bose gas model and 
the PWFRG calculations near the saturation field H 3 for 
some values of the biquadratic interaction f3. (a) /3 = 0, (b) 
P = 1/3, (c) /3 = 0.6. The open circles represent the PWFRG 
results with the retained number of bases m = 60. The solid 
lines show the effective (5-bose gas model with the coupling 
c = —(3/3 + l)//3. In (b) and (c), we draw the free fermion 
curves, which correspond to |c| = oo, as the broken lines for 
comparison. 

FIG. 2. Comparisons of the effective (5-bose gas model and 
the PWFRG calculations for various S at /3 = 0. The dia- 
monds, squares, triangles and circles represent the PWFRG 
results for S = 1/2, 1, 3/2 and 2 respectively. The solid lines 
represent the effective (5-bose gas model with c = 2/(S — 1). 



FIG. 3. Comparisons of the effective 5-bose gas model 
and the PWFRG calculations near the lower critical field 
H c (= A), (a) The pure AF Heisenberg point /3 = 0. (b) The 
AKLT point p = 1/3. The open circles represent the PWFRG 
results with the retained number of bases m = 100. The 
solid lines show the least-square fitting results of the form: 
H = A + qtt 2 M 2 (1 + 7 M + SM 2 ). 



FIG. 4. The H - M curve at (3 = 0.4(« (3 C ) near the 
lower critical field H c (= A). The open circles represent the 
PWFRG result with m = 100. The solid line shows the 
least-square fitting result of the form: H = A + AM 4, + BM 5 . 

FIG. 5. The H - M curve at /3 = 0.6 near the lower crit- 
ical field H c (= A). The open circles represent the PWFRG 
result with m = 100. The solid line shows the least-square 
fitting result. We see the bose gas coupling constant takes the 
negative value. 
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